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ABSTRACT 

We build a model for the density and integrated Sachs- Wolfe (ISW) profile of supervoid and su- 
percluster structures. Our model assumes that fluctuations evolve linearly from an initial Gaussian 
random field. We find these assumptions capable of describing N-body simulations and simulated 
ISW maps remarkably well on large scales. We construct an ISW map based on locations of super- 
structures identified previously in the SDSS Luminous Red Galaxy sample. A matched filter analysis 
of the cosmic microwave background confirms a signal at the 3.2 — cr confidence level and estimates 
the radius of the underlying structures to be 55 ± 28h^^Mpc. The amplitude of the signal, however, 
is 2 — cr higher than ACDM predictions. 

Subject headings: large-scale structure of universe — cosmic microwave background — methods: sta- 
tistical 



1. INTRODUCTION 

Since the first sky maps of the Wilkinson Microwave 
Anisotropy Probe (WMAP) were published, there have 
been claims for the existence of circular features (spots 
and rings) in the c osmic microwave back ground (CMB) 
(|Cruz et al\ [20051: IGranett et all l2008aD . Their origin 
and statistical significance are still debated. Sources, 
such as foreground contami nation, integrated Sachs- 
Wolfe (ISachs and Wolfd[T967t ISW) effect and the more 
exotic conformal cyclic cosmology or cosmic texture have 
been considered viable candidates to explain circular 
CMB anomalies (IRudnick a/.l 120071: IInoue etaZ.1 120101: 
iCruz et a/]|2007t iGurzadvan and Pen rose 201^. A no- 
table feature is the cold spot ^Crva et a/.n200 a). which 
has a mean temperature of —70fj,K in a 5°-radius aper- 
ture. Additionally, lO^K hot and cold spots have been 
identified on 4° scale s associated with super structures 
(jCranett et a/.ll2008aD . 

In this paper we focus on the integrated Sachs- Wolfe 
effect and investigate whether it is possible that large 
spherical fluctuations (supervoids and superclusters) in 
the dark matter field produce the aforementioned fea- 
tures in the CMB. The large-scale ISW effect is expected 
in a universe with accelerated cosmic expansion arising 
either from a dark energy term in flat cosmological 
models, or from spatial curvature. The effect is sensitive 
to both the expansion history and the rate of structure 
formation and provides cons traints on alternative cosmo- 
logical models (e.g. Giannantonio et a/.,,2 008bl as well as 
initial non-Gaussianitv fe.g. lAfshordi and ToUev 2008). 
Cross correlating a galaxy catalog and a CMB temper- 
ature map is the standard way of stu dying the ISW 
signal and it has an extensive literature (jScranton et al\ 



200§, 'Afshordi et aV '2004'; 'Padmanabh an et al\ [2005 : 
Raccanelli et at. ,2008: , Giannantonio 'et al\ 12008a : 
Ho et al.\ l2008riSawangwit et al.\ \20]W ~ Our work is 



similar considering what we measure is related to the 
correlation between the dark matter and the CMB. 
However, we only focus on parts of the sky where the 
signal is expected to be large, regions corres ponding to 
supervoids and supercluster s. Since reports (jCruz et al\ 
120051: IGranett etM l2008aD indicate that the scale of 
these regions is beyond nonlinear scales, our model for 
their average density proflle and ISW imprint is derived 
from the statistics of the linearly evolving prim ordial 
Gaussian density he ld (IPapai and Szapudill2010D . This 
is in contrast with llnoue et al\ ( 2010[ ). who assumed 
a top hat density profile, which is the asymptotic 
final state of a void with steep initial density profile 
(iSheth and van de Wevgaertll20Q4[ ). 

We use the publicly available Hubble Volume Sim- 
ulation of the V irgo Supercomputing Consortium 
(jColberg et a^.l 12000: HVS) to test the Gaussian model. 
This is the N-body simulation with the largest volume, 
which is relevant considering that the gravitational po- 
tential is correlated more strongly than the density. We 
simulate ISW maps by ray tracing through the potential 
and calculating the linear part of the ISW effect. These 
are compared to partial ISW maps generated from sets 
of spherical regions based on our model. 

After gaining some confidence by studying simula- 
tions we create an ISW map fr om real data. We se - 
lect locations on the sky based on lGranett et ali (|2008a'). 
They compiled a list of supervoids and superclusters 
found in the Sloan Digital Sky Survey (SDSS) Luminous 
Red Galaxy (LRG) s ample. The list can be found in 
IGranett et ai] (|2008b[ ). We build an ISW map by plac- 
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ing the theoretical profiles to the given R.A., decl. co- 
ordinates. This map is fitted to a WMAP temperature 
map. 

The structure of the paper is the following: in Sec- 
tion [2] we measure the expected density and ISW profile 
of spherical dark matter fiuctuations in N-body simula- 
tions; in Section [3] we apply the matched filter technique 
to detect the signature of superstructures in the CMB; in 
Section|4]we discuss our results and views of the relation- 
ship between the linear ISW effect and circular features 
on the CMB. 

2. PROFILES IN SIMULATIONS 

In this section we calculate the full expected density 
profiles and ISW imprin ts of spherical overd e nsitie s and 
underdensities based on iPapai and Szapudil (|2010( ) and 
compare them with measurements from the HVS. 



2.1. 



Density Profiles of Superstructures in N-body 
Simulations 



According to IPapai a nd Szapudi (2010), if the aver- 
age density contrast in a Gaussian random field inside a 
sphere of radius R is given, we expect the density con- 
trast at radius r from its center to be: 



(1) 



where (...) stands for ensemble averaging, S{r) is the den- 
sity contrast at radius r and 6in{R) i s the average of 5 
inside radius R. In the Appendix of iDekell ()1981f ) the 
author derived a similar formula for a single, point-like 
location with given density, which is the i? — > limit. 
Equation ([T]) is essentially the two-point function multi- 
plied by a normalization constant. A more general cal- 
culation for the density profile around loc al maxima in a 
Gauss ian random field was carried out by iBardeen et al\ 
(fT986l) 

In order to measure ((5(r))^ one needs to select 
regions with 5in [R) and calculate the average over these. 
However, according to Equation ([T]), ^^^^'j^.-^ gives an un- 

{s(r)5UFt)) 



biased estimator of ■ 



the shape of the profile: 
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We 
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is independent of dm- The 



where we shorten 6{r) and 5in{R) as 5 and (5; 
also use the fact that by definition ((5(r)) 

argument above shows that it is needless to search for 
spheres with a certain Sin{R) to verify Equation ([T]). In 
the rest of this Subsection we deal with the integral of 
d{r) instead of S{r), because the former is measurable 
directly: 
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Fig. 1. — The integral of the spherically averaged density profile 
measured in the HVS and 100 2LPT simulations. Linear theory 
given by Equation l(2)l appears to be a good approximation. 

The linear size of the HVS is 3000h-iMpc. We 
chose to place spheres at 25^ evenly spaced coordi- 
nates. We counted galaxies in concentric spheres with 
radn of lOh-^Mpc, 20h~iMpc, 160h~iMpc and sub- 
tracted the average to get M(r). We normalized this 



as 



M{r) 47r 
M{R.} 3 



M(r) = 47r drr^d{r). 



(3) 



R'^ and calculated the average of the 25^ in- 
stances to get an estimate of the integral of Equation 
As R changes, technically, by weighting M{r) differ- 
ently, one can get the profile of linear superstructures of 
different sizes. A nuisance, we have to deal with, is that 
the estimator given by Equation ([2]) is unstable because 
Sin can be arbitrarily small. The ratio of two Gaussian 
variables with zero means has a Lorentzian distribution, 
which has an undefined variance. This problem can be 
avoided if cases, when \din\ is under a threshold, are ig- 
nored. This constraint leaves the estimator unchanged as 
it can be seen from Equation ([2|). We chose the thresh- 
old based on the variance of Sin in the particular cat- 
alog, e. q. 2^ i^Sfn)- The result is not sensitive to 

the exact choice but the right choice can decrease the 
variance of the estimator substantially. In Figure [T] we 
plotted the measured profiles and their theoretical coun- 
terparts for R — 60h~^Mpc on the upper panel and for 
R = lOOh^^Mpc on the lower panel. We obtained the 
error bars by repeating the measurement on one hundred 
mock catalogs generated with a se cond order Lagrangian 
(2LPT) code (jCrocce et a/.ll2006[ ). The mocks were set 
up to have the same cosmological parameters, size and 
density as the HV simulation. To calculate the theo- 
retical profiles of Equation ([T]) we used a matter power 
spectrum calculated via CAMB (iLewis a/.ll2000( ). The 
measured profiles, both in 2LPT and the HVS, appear to 
be in good agreement with theory within the uncertainty 
up to 400h~^Mpc, the largest radius in our measurement. 

2.2. The ISW imprints of superstructures 

In IPapai and Szapudil (|2010D the ISW effect in the di- 
rection of the center of a superstructure was calculated. 
This was done simply by calculating the spherically sym- 
metric potential of the superstructure, then computing 
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the following integral a long the path of the CMB photon 



tne lollowmg integral a le 
(|Sachs and Wolfeill967D : 



$(r) 



AT 
T 



. aj'(r(r),r) 
Or ' 



(4) 
(5) 



where r denotes conformal time. Similarly, this calcula- 
tion can be carried out for any direction easily in order to 
predict the full ISW profile. In linear theory the potential 
at a a particular redshift is just $(r, z) = $(r, z — 0)-^^. 

We leave the study of nonlinear corrections for future 
work. Since these integrals are linear, the linear ISW 
signal is predicted to be proportional to the average den- 
sity of the superstructure, Sm ■ In this paper we place the 
centers of the superstructures aX z = 0.52, which is the 
median redshift of the galaxy catalog we use in Section 

El 

2.3. Fluctuations in the Potential on the Largest Scales 
and the ISW Profile 

To test the ISW profile we traced rays through the 
HV simulation along the z-axis. We calculated S on an 
800^ grid and calculated $ at z = by using the simple 
formula: 



(6) 



Since we aimed to analyze imprints of superstructures 
at redshift 0.52, we defined the boundaries of the inte- 
gral of Equation ([5]) in a way that put the 25'^ prese- 
lected locations (see Subsection l2.ip at this redshift. We 
used periodic boundary conditions and we integrated up 
to SOOOh^^Mpc from starting points ensuring that 25^ 
of the locations were at redshift 0.52 each time. This 
yielded 25 ISW images of the HVS in the x — y plane. 
Each of these images contains imprints of 25^ superstruc- 
tures. After this, one can follow the procedure discussed 
for density profiles in Subsection 12.11 If the ISW profile 
of a superstructure is normalized by the average density 
of the superstructure, the result will be an estimator of 
the shape of the ISW profile: 



(7) 



The arguments are r2D, the radius in the x — y plane and 
R, the radius in which we know the galaxy count in the 
corresponding volume. 

Here, we would like to remind the reader that 
refers to ensemble averaging. For density, averaging over 
a relatively large volume proved to be a sufficient substi- 
tute in Subsection 12.11 Despite the fact that the cosmic 
variance is large on the largest scales, these modes con- 
tribute only a little to the density distribution since their 
amplitudes are small. (The power spectrum goes to zero 
with k.) For the potential, this picture changes due to 
the 1/fc^ factor. As a result, the ISW effect has a large 
cosmic variance regardless of the large size of the simu- 
lation. To demonstrate this, we projected the density to 
the x — y plane. We split this projection in the middle 
of the X-axis to create two 1500 X 3000Mpc^h~^ areas. 
In Figure [5] in the top row we plotted the histograms 
of these. We did the same with one of the 25 ISW im- 
ages, which is essentially the projection of the potential. 
Those histograms are shown in the middle row. The 
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Fig. 2. — Each row shows histograms of a certain projection 
of density along the z-axis. A simple integral of f 5{x, y, z)dz is 
presented in the first row, an ISW map in the second and its filtered 
version in the third. (See text for details.) The a; — j/ plane is halved 
to test the cosmic variance, hence the two columns. As a result 
of largo fiuctuations in the low fc-modes in the potential, the ISW 
projections show a significant cosmic variance (second row) which 
can be reduced by filtering (third row). 

cosmic variance is significantly larger for the ISW map. 
To demonstrate that the difference is due to modes on 
the largest scales we removed the low fc-modes from the 
potential. After experimenting we found that removing 
modes with |fc| < 27r|;, where L = SOOOh^^Mpc for the 
HVS, gave a visually compelling result (Figure [21 bot- 
tom row). After applying this high-pass filter to a 2- 
dimensional projection, the effect of the cosmic variance 
on the modes with lowest wavenumber will be around 
18%. 

Having removed the low fc-modes from both the density 
and the potential we computed the average in Equation 
([7]). As before, we omitted cases when the denominator 
was too close to zero. (See Subsection 12.11 ) The results 
for R = 60h~^Mpc and R = 100h~^Mpc are shown in 
Figure 131 The error bars show the uncertainty calculated 
from the HVS itself due to the long CPU time required 
to calculate the ISW images. Predictions of the linear 
model are calculated by removing the low fc-modes from 
the power spectrum when doing calculation according to 
Equation ([T|). 

This approach has to be modified slightly when the 
CMB is given and late time anisotropies are affected by 
every mode of the galaxy distribution. The previous ex- 
ercise still proves that low fc-modes are to be ignored if 
we want to study the imprints of superstructures statis- 
tically. This can be achieved by filtering out these modes 
from the CMB before performing such analysis. In ad- 
dition one has to remove the same components from the 
theoretical ISW profile before fitting the data. This is 
shown in Figure IH We filtered the ISW map and we 
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Fig. 3. — The spherically averaged ISW profile from the HVS. 
This is defined by Equation {Tjl. Fourier- modes from the density 
were filtered out above A = 0.2L, where L is the linear size of the 
HVS. The linear approximation is calculated as described in the 
text. See Equations Q — ISJ. 
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Fig. 4. — The same as Figure|3]with the exception that the ISW 
map was filtered, not the density. The error bars represent the 
uncertainty of the measurement. 

filtered the ISW profile, which is implicitly expressed as 
f{'r2DiR), on the right side of Equation ([T]), so that only 
modes with \k\ > 2ttj^ remained. After these manipu- 
lations the measurement remained consistent with linear 
theory. 

The success of t hese tests proves the vali dity of the 
Gaussian model of iPapai and Szapudil (|20100 . We stud- 
ied Euclidean ISW projections because or their simplic- 
ity but there is no reason to think that our findings are 
restricted to the Euclidean case. The only major differ- 
ence is that on a sphere one deals with spherical harmon- 
ics and instead of low A;-modes we speak of low Z-modes 
analogously. 

3. ISW MAP OF SDSS SUPERSTRUCTURES 



Having calculated the full linear ISW profile of super- 
structures in the dark matter density, we attempt to de- 
tect their presence on the CMB. With the locations of 50 
superclusters an d 50 supervoid s found in the SDSS DR6 
LRG sample by 'Grane tt et oD ()2008al ) we created ISW 
maps and used the matched filter technique to measure 
the amplitude of the ISW signal. Our ISW maps are 
statistical in the sense that the individual ISW profiles 
were predicted by using statistical properties of Gaus- 
sian, spherically symmetric density fiuctuations rather 
than the particular (poorly observed) realization of the 
density field. 

3.1. Superclusters and Supervoids in the SDSS LRG 
sample on 100h~^Mpc scales 

Large regions with significantly high overdensities or 
low underdensities are called here superclusters or su- 
pervoids. In a given a dataset these regions can be de- 
fined in many diffe rent way s . Void finder algorithms were 
compared in IColberg et a/.l ()2008f ) . Each implementation 
is based on a certain notion of a void. Voids can have 
spherical or irregular shapes, finders can be parameter 
free or there can be a hard coded density threshold or 
radius to describe what a void is like. It is important 
to emphasize that despite the similarities, any analysis 
using the output of a void finder will be sensitive to its 
particulars. This renders our study qualitative in nature. 

Our work is based on supervoids and superclusters 
found in a subsample of the SDSS DR6 LRG catalog by 
ZOBOV and VO BOZ, p arameter free v oid and cluster 
finder algorithms (jNevri nck 2008; Ncvrin ck eraLl l2005'). 
ZOBOV and VOBOZ calculate the Voronoi tessellation 
of the data in order to approximate the density field with 
a simple function, a funcion that is constant inside poly- 
hedra. Voids are found by using the concept of drainage 
basins from Geography. Clusters are identified analo- 
gously. 

A table with the propert ies of 50 supervoids an d 50 su- 
perclusters can be found in lGranett et al\ ()2008b[ ). Their 
sizes and densities cannot be quantified simply due to 
their shapes and survey boundary effects. The bias, red- 
shift error and the shot noise are also sources of uncer- 
tainty for both size and density. This is why we used only 
the directions of the centers of these superstructures in 
creating ISW maps in Section [3?3l The error in the an- 
gular coordinates is relatively small. 

3.2. CMB and ISW Maps 

We used the Internal Linear C ombination (ILC) Map 
from the WMAP 7-year release ()Gold et aLl 120101 ) with 
the KQ75y7 galactic foreground and point source mask. 
We employ ed Hcalpix pixelization at resolution Uside = 
64 (|G6rski et al. .2005i) . In general, we strove to carry 
out the measurement with up t o date maps and ma sks 
while remaining consistent with iGranett et al\ ()2009fl . 

Section[2]concluded that a Gaussian model gave a good 
approximation to the expected ISW signal from a large 
spherical region in the dark matter field with given den- 
sity contrast. It also became apparent that in an anal- 
ysis of statistical ISW maps low Z-modes were to be fil- 
tered out due to the large cosmic variance of the poten- 
tial affecting these modes. This also lowers the cosmic 
variance. The lowest mode to be taken into account is 
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minimum of Equation ([5]) : 





n,n = 
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Fig. 5. — The expected ISW imprint of a spherical region with 
radius of 120h~^Mpc and an average 5 of 1 at redshift 0.52. Each 
curve shows the same profile without modes under their respective 



determined by the sky coverage of the ISW map. Experi- 
menting with the HV Simulation gave us kmin ~ 27r-^ for 
an Euchdean ISW map, where L is the hnear size of the 



map. This translates into Iminih: 



1 



4-7r 
A 



, where A 



is the angular size of ISW map. The SDSS DR6 imaging 
area about the north galactic pole covers roughly 8000 
square degrees which gives us Imin = 12. 

Due to the large uncertainty in their size and density 
we assume that each supervoid is a realization of the 
same statistical entity when we creating ISW maps. We 
consider superclusters similar in nature but with an op- 
posite density contrast. 

An ISW map is made simply by placing the expected 
ISW profiles o f supers tructures to the locations given in 
iGranett et al\ (|2008bf ). We set Sin to 1 for clusters and 
— 1 for voids in formulae of Section [2] For every radius 
there is a corresponding ISW map. We calculated the 
aim coefficients of the spherical expansion of both the 
ISW maps and the CMB map, and zeroed the modes be- 
low Imin — 12 before transforming the a/„i's back. These 
manipulations we re carried out with routines from the 
Healpix package (iGorski et al\ 120051 ). Figure [5] shows 
how the ISW profile changes when modes are removed. 

3.3. The Matched Filter Technique 

We continue with measuring the amplitude of the ISW 
maps created in Section 13.21 in the WMAP ILC map. 
To do this we used t he same matched filter technique as 
IGranett et al\ ()2009[ ). A likelihood function is defined 
as: 

L(X) = —-{Tcmb — XTlsw)C~^{Tcmb — AT/sw) , (8) 

where Tcmb is a foreground cleaned CMB map, in this 
case the WMAP ILC map, and Tisw refers to the ISW 
maps. The pixel covariance matrix, C, is computed from 
the best fitting ACDM power spectrum. This likeli- 
hood is based on the assumption that the Tjsw and the 
primary CMB are the realizations of two very different 
Gaussian random fields and the latter is given by the 
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(9) 
(10) 

As explained in Section [3.2[ we filtered out modes up to 
Imin = 12 from both Tjsw and Tcmb, and before building 
the covariance matrix from C/'s we set Cq = Ci = ... = 
Cii=0. 

Since the exact inverse of the covariance matrix is un- 
stable for partial sky maps, we regularized C by calcu- 
lating a pseudo-inverse omitting the noisiest modes. The 
likelihood function above defines an error on A* as 



AA2 



(11) 



3.4. Results 

We calculated A* according to Equation Since we 
have ISW maps for different void/cluster radii, the result 
is a function, A(i?). This is presented in Figure [6] in the 
bottom row on the left panel. The error bars are given 
by Equation (|11|) . From now on every time we refer to 
a A we mean the value at the minimum of the likelihood 
function, so we drop the superscript *. If we marginalize 
\{R) / AX{R) over R we get a 3.24-(t detection. Similarly, 
the mare inalized value of i? is 55 ± 28h~iMpc. 

A(i?) can be interpreted as the average density contrast 
in one hundred superstructures inside radius i?, where 
the supervoids are taken into account with a minus sign: 

(12) 

where superscripts 50c and 50v refer to the average of 
50 superclusters and 50 supervoids. This interpretation 
originates from the linear relationship between the ISW 
signal and the density fluctuations, and our choice of 
Sin = ±1 when building ISW maps in Section [221 In ad- 
dition to the data curves, in Figure [S] we plotted the pre- 
diction of linear theory for the average of the 50 largest 
superclusters. We assumed a Gaussian Probability Dis- 
tribution Functio n (PDF) which is reason able in light 
of the findings of P"apa i and SzapudU ()2010( ) . The vari- 
ance of the PDF was calculated by CAMB ([Lewis et al\ 
[2000h using the late st WMAP cosmological parameters 
( Jarosik et al\\20Wt l and it was scaled with the growth 
function to the median redshift of the superstructures 
(z = 0.52). The number of independent clusters was es- 
timated as Vsurvey/Vcluster whcrC Vsurvev I S the V Olume 

of the LRG sample used by IGranett et al\ ([2009D . A is 
consistently several times larger than the predicted Sin- 
As a consistency test we split the voids and clusters into 
two groups. Group 1 consists of the first 25 voids and 
25 clusters with the highest significance in the supple- 
mentary tables of IGranett et al\ (|2008bl) . The other 25 
voids and 25 clusters make up group 2. We created two 
ISW maps, one from each of the groups. After repeating 
the whole procedure for these, we got the center and the 
right panel in the bottom row of Figure [6] As expected, 
superstructures associated with lower significance pro- 
duce lower signal on the CMB and allow larger deviation 
from the average A. This supports the hypothesis that 
the signal is due to the ISW effect. 

We also checked the sensitivity of the results to the 
choice of Imin, the lowest spherical harmonics included 
in the analysis. The result for different choices of Imin 
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Fig. 6. — The top and the mid row contain contour plots of the best fitting amplitude of the ISW maps, A, and their relative uncertainty, 
<T, (see text for their definition) as functions of Imin and R. The bottom row shows cross sections of the contour plots from the top row at 
Imin = 12, as well as the prediction for the average density contrast of superstructures of a certain size. In the left column the fitting is 
for an ISW map consisting of the 50 most significant supervoids and 50 superclusters. The middle column shows the same for best 25-25 
while the right for the remaining 50 superstructures not included in the middle column. 



is shown in the mid and top rows of Figure [51 The 
top row shows X{R,lmin), while the mid row shows 

0-{R,lmin) = X{R, Imin) /^\{R, Imin) On COntOUr plotS. 

The resuh seems robust for wide range of Imin^. 

We performed another consistency check to acquire a 
better understanding of the uncertainty of A and its re- 
lationship with 6in- We created ISW maps from single 
superstructures. We used the matched filter technique 
to determine a A and its error as before. The average A 
for voids and for clusters separately are plotted in Fig- 
ure [71 In this case we defined the uncertainty of A as 
the standard deviation of its mean. We also plotted the 
previous results for the previous combined ISW map. All 
of this is for Imin = 12. The three measured curves are 
consistent with each other, the error bars show the two 
a deviation. The single void and cluster fit allow wider 
range of A than the combined fit. 

4. DISCUSSION 



This paper is the conti nuation of the work started in 
iPapai a nd Szaptidl (|20100 . Our goal is to estimate the 
ISW imprint of large spherical dark matter overdensities 
and underdensities. In Section[21 first, wc tested a simple 
Gaussian model for density profiles in the HVS. We cre- 
ated an ISW map by ray tracing through the simulation 
and computing the linear ISW effect. We neglected the 
nonlinear part of the ISW effect, as it is small in compar- 
ison at z = (jCai et aLI 12010 V Another simplification 
was that our ISW map was Euclidean. For the purpose 
of testing this has no relevance. 

As it can be seen from Figure [I] the density profiles 
followed linear theory given by Equation ^ within l-cr 
uncertainty up to 400h~^Mpc, the largest scale in our 
study. This is not surprising, since the density profile is 
equivalent of the two-point function and on large scales 
high-order clustering is not important. When we aver- 
aged the density profiles we only used locations where 
the average density contrast in a certain radius was in the 
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Fig. 7. — After fitting simple ISW maps made from the expected 
imprint of single superclusters or supervoids we were left with SO- 
SO best amplitudes, A's. Their averages with the 2(t uncertainty 
are plotted with the amplitude of a full ISW map consisting of all 
100 imprints. These measurements were done with Imin = 12. The 
theoretical prediction (continuous line) for the density contrast is 
plotted along. 

2-(7 (positive or negative) tail of the PDF. This demon- 
strates that not even the extreme cases are affected by 
nonhnearities significantly. 

We demonstrated that cosmic variance was a much 
more important factor for the ISW effect than for the 
density. Since cosmic variance is large for low fc-modes 
and the potential is non- vanishing as ^ 0, these modes 
have relatively large effect on the ISW profile. Surely 
large variance modes should be removed from the CMB 
and the ISW maps if the ISW maps are based on the 
statistics of density fluctuations and not on the observed 
density. First we measured the ISW profiles after filtering 
out the low fc-modes from the HVS density and conse- 
quently the ISW map. Filtering the two-point function 
in the same way gave a good match to the data. (See 
Figured) 

If the ISW map is given, as in the CMB, another route 
has to be taken. We filtered the complete ISW map and 
compared it to the theoretical profile after removing the 
same modes. The result is shown in Figure U 

In Section [3] we used matched filter technique to de- 
tect the ISW signal of superstructures in the CMB. We 
closely followed the procedure described in lGranett et al\ 
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((20091) . The major difference is that our ISW map was 
built from theoretical assumptions about the shape of 
ISW profiles, while in iGranett et al\ ()2009l ) the authors 
used an analytic curve which fitted the measured pro- 
files the best. Because of this, the significance of our 
measurement can readily be interpreted. We estimated 
the marginalized significance of our measurement to be 
(7 = 3.24. (See Subsection [3]4] for details.) The interpre- 
tation of the best fitting amplitude of the ISW map is still 
not without difficulties. It is several times higher than 
the anticipated signal (Figures [H] and [7]), although it ap- 
pears to be 2-(T higher than ACDM predictions. Despite 
the difference in amplitude, the theoretical and the mea- 
sured curves run parallel with each other, which supports 
that the signal is related to the ISW effect. We caution 
though that the uncertainty can still be underestimated, 
since it was derived from data alone. 

The ISW effect provides an explanation for another 
feature of observations. On average supervoids tend to 
have a hot ring around an i nitial cold dip and th e oppo- 
site is true for superclusters (jOranett et aLl[20Q9l ). When 
we erased large scale fluctuations in the microwave back- 
ground, the expected ISW profile changed sign as seen 
in Figure [5] This means that in certain, not particularly 
unique realizations the average profiles of 50 supervoids 
or superstructures from a finite area of the sky can have 
rings in the ISW context, simply due to the cosmic vari- 
ance of low fc-modes. 

Overall we can say that the shape of the measured sig- 
nal follows the predicted ISW profile while its amplitude 
exceeds expectation. This is a good reason to investi- 
gate further by studying galaxy surveys other than the 
SDSS. A key to a quantitative study is the well-measured 
galaxy density. Since the model described in this paper 
is a model for density fluctuations around a particular lo- 
cation, which is not necessarily a maximum, it is enough 
to know the density at this location accurately. 

We thank Mark Neyrinck for sharing his knowledge on 
void and cluster finders. The authors were supported by 
NASA grants NNXI0AD53G and NNG06GE71G, and 
the Polanyi program of the Hungarian National Office 
for the Research and Technology (NKTH). 

The Hubble Volume Simulations were carried out by 
the Virgo Supercomputing Consortium using computers 
based at the Computing Centre of the Max-Planck So- 
ciety in Garching and at the Edinburgh parallel Com- 
puting Centre. The data are publicly available at 
|http: / /www.mpa-garching.mpg.de/NumCos[ 



Crocce, M., S. Pueblas, and R. Scoccimarro, 2006, MNRAS 373, 
369. 

Cruz, M., E. Martinez-Gonzalez, P. Vielva, and L. Cayon, 2005, 

MNRAS 356, 29. 
Cruz, M., N. Turok, P. Vielva, E. Martinez-Gonzalez, and 

M. Hobson, 2007, Science 318, 1612. 
Dekel, A., 1981, A&A 101, 79. 

Giannantonio, T., R. Scranton, R. G. Crittenden, R. C. Nichol, 
S. P. Boughn, A. D. Myers, and G. T. Richards, 2008a, 
Phys. Rev. D 77(12), 123S20. 

Giannantonio, T., Y. Song, and K. Koyama, 2008b, Phys. Rev. D 
78(4), 044017. 

Gold, B., N. Odegard, J. L. Weiland, R. S. Hill, A. Kogut, C. L. 
Bennett, G. Hinshaw, X. Chen, J. Dunkley, M. Halpern, 
N. Jarosik, E. Komatsu, et al, 2010, ArXiv e-prints 1001. 45S5. 



8 



Gorski, K. M., E. Hivon, A. J. Banday, B. D. Wandclt, F. K. 

Hansen, M. Rcincckc, and M. Bartelmann, 2005, ApJ 622, 759. 
Granett, B. R., M. C. Ncyrinck, and I. Szapudi, 2008a, ApJ 683, 

L99. 

Granett, B. R., M. C. Ncyrinck, and I. Szapudi, 2008b, ArXiv 

e-prints 0805.2974. 
Granett, B. R., M. C. Ncyrinck, and I. Szapudi, 2009, ApJ 701, 

414. 

Gurzadyan, V. G., and R. Penrose, 2010, ArXiv e-prints 
1011.3706. 

Ho, S., C. Hirata, N. Padmanabhan, U. SeljaJc, and N. Bahcall, 
2008, Phys. Rev. D 78(4), 043519. 

Inoue, K. T., N. Sakai, and K. Tomita, 2010, ApJ 724, 12. 

Jarosik, N., C. L. Bennett, J. Dunkley, B. Gold, M. R. Greason, 
M. Halpern, R. S. Hill, G. Hinshaw, A. Kogut, E. Komatsu, 
D. Larson, M. Limon, et al, 2010, ArXiv e-prints 1001.4744. 

Lewis, A., A. Ciiallinor, and A. Lasenby, 2000, ApJ 538, 473. 

Ncyrinck, M. C., 2008, MNRAS 386, 2101. 



Ncyrinck, M. C, N. Y. Gncdin, and A. J. S. Hamilton, 2005, 

MNRAS 356, 1222. 
Padmanabhan, N., C. M. Hirata, U. Seljak, D. J. Sclilegel, 

J. Brinkmann, and D. P. Schneider, 2005, Phys. Rev. D 72(4), 

043525. 

Papal, P., and I. Szapudi, 2010, ApJ 725, 2078. 

Raccanelli, A., A. Bonaldi, M. Negrello, S. Matarrese, G. Tormen, 
and G. de Zotti, 2008, MNRAS 386, 2161. 

Rudnick, L., S. Brown, and L. R. Williams, 2007, ApJ 671, 40. 

Sachs, R. K., and A. M. Wolfe, 1967, ApJ 147, 73. 

Sawangwit, U., T. Shanks, R. D. Cannon, S. M. Groom, N. P. 
Ross, and D. A. Wake, 2010, MNRAS 402, 2228. 

Scranton, R., A. J. Connolly, R. G. Nichol, A. Stebbins, 
L Szapudi, D. J. Eisenstein, N. Afshordi, T. Budavari, 
I. Csabai, J. A. Prieman, J. E. Gunn, D. Johnson, et al., 2003, 
ArXiv Astrophysics e-prints arXiv:astro-ph/0307335. 

Sheth, R. K., and R. van de Weygaert, 2004, MNRAS 350, 517. 



